Tricritical point in heterogeneous fc-core percolation 
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/c-core percolation is an extension of the concept of classical percolation and is particularly relevant 
to understand the resilience of complex networks under random damage. A new analytical formalism 
has been recently proposed to deal with heterogeneous /c-cores, where each vertex is assigned a local 
threshold ki. In this paper we identify a binary mixture of heterogeneous /c-cores which exhibits a 
tricritical point. We investigate the new scaling scenario and calculate the relevant critical exponents, 
by analytical and computational methods, for Erdos-Renyi networks and 2d square lattices. 



Tricritical points (TCPs) are interesting critical phe- 
nomena in statistical physics, constituting a natural 
switch between 1 st and 2 nd order phase transitions. In 
other words, a TCP affords the possibility to smoothly 
control the order of the phase transition by tuning the 
appropriate parameter. In the language of critical phe- 
nomena, this is equivalent to identifying the position of 
a phase transition from the extrapolation of the order 
parameter- while usually impossible in 1 st order tran- 
sitions, it can be done in the case of continuous ones 
PQ. In the case of percolation models, there has recently 
been an attempt to identify a TCP in a model mixing 
elements of classical and explosive percolation in a lattice 
[2], although there is now evidence that the explosive per- 
colation transition is continuous |3J. Other recent models 
which allow control of the order of the transition include 
explosive percolation on scale free networks [4], and de- 
pendency groups on interdependent networks [5]. In this 
paper we establish, for the first time, the presence of a 
TCP in a simple extension of classical percolation, namely 
heterogeneous k-core (HKC) percolation, which has the 
advantage of a sound analytical approach on random and 
complex networks [6j [7] . We show analytical evidence of 
a TCP in Erdos-Renyi graphs and numerically we find 
similar phase diagram topology in the square lattice. Fi- 
nally, our model appears to be in the same universality 
class of a model which reproduces non-trivial signatures 
of liquid-glass transitions, including the higher-order glass 
singularity predicted by mode-coupling theory |8J. 

A k-coie is defined as the maximal network subset 
which survives after a culling process which recursively 
removes all the vertices (and adjacent edges) with less 
than k neighbors. As a generalization of the concept of the 
giant component, the k-core gives a deeper insight into the 
structure and organization of complex networks. It has 
been thoroughly investigated on Bethe lattices |9J, random 
graphs [6j[T0] and, using a numerical approach, on various 
types of lattices [llj. The k-core percolation analysis 
has found several applications in varied areas of science 
including protein interaction networks [12], jamming [T3] , 
neural networks [14], granular gases [15] and evolution 
[T6] . Important insights into the resilience of networks 



under damage [17] and spreading of influence in social 
networks [18J are gleaned from an understanding of the 
k-core structure of the network. As in Ref. [SI HI], we 
can study k-core percolation on networks after randomly 
removing a fraction 1 — p of vertices. We use the treelike 
properties of the configuration model [20J, in which the 
number of loops vanishes as N — » oo, which guarantees 
that if a k-core exists, it must be infinite, at least if 
k ^ 2 [6[ [9] . In the HKC extension [7] each vertex has its 
own threshold and the culling process is based on local, 
vertex-dependent rules. Although Baxter et al. developed 
results for an arbitrary distribution of vertex thresholds, 
they study binary mixtures of vertices of types a and 6, 
with thresholds k a = 1, ^ 3. The first heterogeneous 
models of this kind were investigated by Branco [21 J on a 
Bethe lattice, whereas the related problem of bootstrap 
percolation (BP) has been much studied on regular lattices 
[TT | [22 ] . Here we focus on the case k = (fc a , fc&) = (2, 3). 



We start with a binary mixture (k a ,kb), where vertices 
have been randomly assigned two thresholds k a and 
(say k a < kb) with probability r and 1 — r, respectively. 
Finite clusters are a possibility when k a = 1 and so we 
must make a distinction between M a ^, the probability that 
a randomly chosen vertex belongs to the HKC, and *S a &, 
the probability that a randomly chosen vertex belongs 
to the giant component of the HKC. We will show that 
in the case k = (2, 3) these two quantities are coincident, 
but there are relevant examples where they are not [7J. 



In the original k-core formalism, given the end of an 
edge, a (k — l)-ary subtree is defined as the tree where, as 
we traverse it, each vertex has at least k—1 outgoing edges, 
apart from the one we came in. Instead, considering a 
HKC, every vertex i may have a different threshold ki. 
The (k{ — l)-ary subtree, then, is the tree in which, as we 
traverse it, each encountered vertex has at least k{ — 1 child 
edges. We define Z as the probability that a randomly 
chosen vertex is the root of a (ki — l)-ary subtree. Taking 
advantage of the local treelike nature of the configuration 
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model, Z is related to as [7]: 

oo 

M ab (p) = M„(p) + M 6 (p) = pr ^ P(g)$J»(Z,Z)+ 

q=k a 

oo 

+p(l-r)2^(9)*J i (^^) (!) 

where M a (6)(p) is the fraction of nodes of type a(b) in the 
HKC, respectively, P(q) is the degree distribution and we 
have used the convenient auxiliary function: 
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The quantity &q a(b) (Z, Z) in represents the probability 
that a vertex of type a(b) of degree has at least k a (kb) 
edges which are roots of a (ki — l)-ary subtree. This quan- 
tity is summed over all possible degrees, taking account 
of the relevant vertex type fraction. The self-consistent 
equation for Z is: 

qP(q). 



Z=pr J2 



+p(i 



(q) 



q=k b 



qP(q) 
(q) 



^(ZZ) 



(2) 



We now consider the probability X, that a randomly 
chosen edge leads to a vertex which is the root of an 
infinite (ki — l)-ary subtree. In the case of a binary 
mixture, X is written as 17 
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The fraction of vertices in the giant HKC <S a 5, then, is 
given by S a b(p) = S a (p) + Sb(p), where the fraction of 
nodes of type a is S a (p) = prY^Lk* P(q)®q a ( X i z ) and 
an analogous expression holds for Sb(p). 

For k a = 1, kb ^ 3 mixtures on the Bethe lattice, the 
phase diagram shows a critical line which meets a first 
order line at a critical end point and a critical point at the 
end of a two-phase coexistence between a low and a high 
density phase [7]. The two lines do not match at a TCP, 
because the 1-nodes are so robust that a 1-rich phase 
is stable to damage at intermediate compositions even 
when the /c^-rich phase has collapsed. Let us consider now 
the case k = (2, 3), with a degree distribution such that 
Eg Q 2p (q) < oo. We can rewrite Z (Eq. ^ as pf(Z) = 1 
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f(Z) 
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FIG. 1. Phase diagram of the k = (2,3) mixture, showing 
the total mass of the percolating HKC cluster at different 
compositions r, for ER networks with z\ — 10. The TCP at 
r — 1/2 separates a line of 1 st order transitions (dashed) from 
the 2 nd order line (solid). The masses of the 2-rich-core (blue) 
and the 3-rich-core (red) in the giant HKC are also shown. 
The inset shows the phase diagram in the (r,p) space. 



and similarly rewriting Eq. [3] as h(X, Z) = 1/p with 



h(X, Z) = r 
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These two equations differ only in the first (fractionary) 
part of the sum. The X-dependent (positive) general term 
of the series is monotonically decreasing for any < X $C 
1, meaning that Eq. Q has only one non-zero solution 
when Eq. ([2| has a non-zero solution and therefore X = Z 
for the k = (2, 3) mixture (and £23 = M23). We expect 
this property to be true for any mixture with nodes of 
type k ^ 2. 

We now explicitly show that the k = (2,3) mixture 
presents a TCP for an Erdos-Renyi (ER) degree dis- 
tribution with mean degree Z\ P(q) = z\ exp(— Z\)/q\. 
Using the condition X = Z, the equation pf(Z) = 1 
fully solves the problem of finding the onset of the 
giant HKC, and the function f(Z) becomes f(Z) = 



{1 



-z\ Z 



[1 + (1 — r)ziZ]}/Z. It is now clear that 



f'(Z) < for every r > |, implying that the only solution 
is the trivial one Z = 0, with a de-percolating 2 nd order 
phase transition occurring at the critical occupancy prob- 
ability p c = l/rz\. For r < |, f(Z) has a maximum at 
< Zm < 1- This implies the presence of a 1 st order tran- 
sition and a coexistence between a HKC phase of strength 
M23 {Zm), given by ([T]), and the non percolating phase 
at Z = 0. The expansion of f(Z) for r ^ \, Z(p) -» 0, 
as p -> p+, yields f(Z) = rz 1 + (1/2 - r)z\Z + O (Z 2 ) 
showing that the maximum of f(Z) continuously matches 



the Z = line exactly at r t 



where a TCP is present. 
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We show the computed phase diagram of the k = (2,3) 
mixture in Fig. [T] 
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FIG. 2. Re-scaling of the corona mass C23 and the mean 
corona cluster size £23 at the TCP on ER networks. The 
data range in size from N — 2 9 to N — 2 18 via successive 
doublings. We find the exponent ratios f3/v = 0.34(5) and 
sy/v = 0.39(4) from the scaling of C23 and £23 a t Pt (insets) 
and show the data collapse achieved with those exponents 
(main panels N = 2 14 . . . 2 18 ; >, 0, ★,©,□, respectively). 

We now calculate the critical exponents for this mix- 
ture, in particular at the TCP at r t = ^, Pt = j^- The 
expansion of the order parameter M2s(p) for p — » p+ at 
r ^ \ and p —> p* + (the border of the coexistence region) 
at r < \ yields three different values for the exponent /?: 

{2 1/2 < r < 1 
1 r = 1/2 (6) 

1/2 < r < 1/2 

The exponent f3 takes a unique value at the TCP, and 
agrees with the values found by Branco on the Bethe 
lattice [2T| . However, in this work the presence of finite 
size cores had not been properly handled and it was 
erroneously assumed that the phase diagrams of the k = 
(1,3) and the k = (2, 3) mixture had the same topology. 
The exponent = 1/2 for r < 1/2 corresponds to the 
usual hybrid phase transition seen in /c-core percolation, 
a discontinuous transition which combines with critical 
fluctuations (only on the percolating side) as usually found 
in 2 nd order transitions. To our knowledge, the k = (2,3) 
mixture is the first model displaying a TCP adjacent to 
a hybrid phase transition. 

It has been shown that subsets of the HKC called 
corona clusters have the same critical properties of the 
HKC ITU] . The corona vertices have exactly k{ neigh- 
bours in the HKC, and form finite clusters whose mean 
size diverges when approaching the threshold from above. 
The corona clusters provide a more convenient order pa- 
rameter for numerical study of the model on random 
networks, in contrast to the HKC where only one (in- 
finite) cluster survives. Using the configuration model 
with ER degree distribution we simulated the k = (2,3) 
mixture for various sizes. The typical ansatz of finite 
size scaling for a continuous transition is that any quan- 
tity Y scaling as Y ~ (p — p c )~ x should have the form 
Y = N x ^ v F [(p — p c )7V 1/V ] , where v is the correlation 
length exponent and F is some scaling function. Given 
the universal nature of F we expect to see data collapse 



in a plot of YN~ X ^ V against (p - p c )N^ u . Comput- 
ing the mass of the heterogeneous corona 023(h) at the 
TCP for various sizes we find f3/v = 0.34(5) (Fig. [§. 
Similarly for the mean corona cluster size £23, we find 
7/z/ = 0.39(4). We determine the exponent v — 2.86(9) 
by the scaling of the effective percolation threshold with 
size p ave — p c ~ N~ x / V , where we have located p ave 
from the peak of the susceptibility of the corona mass 
AC 23 = ((C 23 ) 2 - (C 23 ) 2 ) 1/2 . We find good data collapse 
with these exponents in the scaling window at the TCP 
(Fig. [2), and fit the exponents /3 = 0.9(90) and 7 = 1.13(1), 
the former being close to the value calculated analytically. 
The behavior of the strength of the HKC along the edges 
of the coexistence region near the TCP for r —> ^ + allows 
us to calculate analytically the subsidiary tricritical expo- 
nent p u defined by M*(r) - (\ - r)^ u [23J. For k = (2, 3) 
we find /3 U = 2. 

The tricritical crossover exponent cp t describes the 
change of the critical line as the TCP is approached [24] . 
Thus, we write the critical line in terms of two scaling 
fields fi± and /in, perpendicular and tangent to the critical 
line, respectively. Given the simplicity of the model, this 
calculation can also be done analytically for ER networks. 
The rotation defining the critical fields is 

( V± \ = ( costf sintf \ (p-pt \ 
y fi\\ J \— sin d cos d J \ r — r t J 

with tan^ = 4/zj.. Close to the TCP, the critical line has 

1 /2 

a behavior ~ ; with a crossover exponent ip t = 2 
(Fig. [T]) . We expect that the above critical behavior (as 
well as the values of the critical exponents) is reproduced 
by all degree distributions with finite second moment. 




FIG. 3. Phase diagram for the k = (2,3) lattice model showing 
the threshold density p against composition r for several sizes 
L = 64(D), 128(o), 256(o), 384(0), 512(x), 768(*), 1024(>). 
The arrows indicate the location of the TCP for each L. The 
inset shows the Binder cumulant U4(p c ) for various sizes, in- 
dicating the narrowing of the tricritical region for increasing 
size. On the right are sample configurations at the threshold 
density in the continuous transition region (top), at the TCP 
(centre) and in the discontinuous region (bottom); in each case 
the critical HKC cluster is red. 
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We simulated the k = (2,3) model on a 2d square 
lattice and located the TCP at small concentrations of 
k = 2 vertex types (Fig. [3]). On the lattice, the analogous 
bootstrap percolation (BP) model has been much studied 
and it is known [25] that for k ^ d + 1 a discontinuous 
transition occurs only at p = 1. For k ^ d the transition 
is continuous, although the critical exponents have values 
which in general depend on k. On the continuous side of 
the TCP we can expect the usual scaling of the threshold 
density p ave —p c ~ whereas for r = we might ex- 

pect the scaling form found in BP p ave —p c ~ 1/ log L [26] . 
although numerical simulations have struggled to confirm 
this scaling in several cases As shown in Fig. [3] the 
TCP moves toward r = with increasing size (determi- 
nation of the precise scaling with L requires far larger 
sizes and is the subject of further work). In fact, there 
is a finite window of r over which the transition slowly 
changes from 1 st to 2 nd order, and this window becomes 
sharper with increasing system size. We quantify this 
with the Binder cumulant U A {p c ) = 1 - (M) 4 /3(M 2 ) 2 
which has the value U4 = | on the 1 st order side and 
on the 2 nd order side (inset of Fig. [3]). Data collapse near 
the TCP does not work due to the presence of different 
scaling regimes. We determined the critical exponents 
at the TCP for the largest size simulated (L — 1024) 
and found the exponents j3 = 0.31(5), 7 = 2.51(3) and 
v = 1.39(9). Exponents 7 and v at the TCP are very 
close to their values for ordinary percolation on a 2d 
lattice (and a little smaller than the ones of explosive per- 
colation |2J). Exponent /?, instead, is significantly larger. 
The fractal dimension of the tricritical HKC clusters is 
D = 2 — f3/v = 1.77(8), somewhat smaller than ordinary 
percolation (D = 1.879), reflecting the presence of large, 
jagged voids in the k = (2, 3) mixture at the TCP. The 
unusual finite-size effects in this model are reflected in a 
violation of the hyperscaling relation. 

In contrast with the k = (2,3) case, the phase diagrams 
for k a = 1, kt> ^ 3 mixtures [7] do not present a TCP. 
Moreover, the analytical properties of f(Z) and h(X, Z) 
indicate that TCPs are also absent in mixtures of type 
k a = 2,kb > 3. Though far from ubiquitous, a TCP 
is indeed present in the k = (2, 3) mixture, not only 
on the Bethe lattice but also in ER graphs and regular 
square lattices. This case appears to be peculiar because 
the resiliences of the 3-rich-phase and the 2-rich-phase 
are sufficiently close that the two phases collapse at the 
same damage fraction, leading to a complete failure of 
the HKC, either through a 1 st or a 2 nd order transition. 
This phenomenon may occur in cases when a mixing of 
a continuously and a discontinuously failing phase is not 
too heterogeneous {k^ = k a + 1). It is intriguing to note 
that the case k = (2, 3) almost exactly maps onto a model 
of glasses, recently studied on the Bethe lattice [8J, which 
appears to be in the same universality class. 

In conclusion, we have presented a new model of HKC 
percolation which supports a smooth interpolation be- 



tween classical percolation and a 1 st order phase transition 
through a TCP. We are able to identify a new tricritical 
scaling scenario and calculate, both by analytical and 
numerical methods, critical exponents which are different 
from the ones of known percolation transitions. We prove 
the presence of this critical phenomenon in ER graphs, 
and we also get strong numerical evidence in the square 
lattice. The capacity to govern the order of phase tran- 
sitions in randomly damaged networks may constitute a 
step towards a more effective infrastructure for network 
protection. 
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